A Necessary and Sufficient Symbolic Condition for the Existence of Incomplete Cholesky Factorization

نویسندگان

  • XIAOGE WANG
  • RANDALL BRAMLEY
  • KYLE A. GALLIVAN
چکیده

This paper presents a suucient condition on sparsity patterns for the existence of the incomplete Cholesky factorization. Given the sparsity pattern P(A) of a matrix A, and a target sparsity pattern P satisfying the condition, incomplete Cholesky factorization successfully completes for all symmetric positive deenite matrices with the same pattern P(A). This condition is also necessary in the sense that for a given P(A) and target pattern P, if P does not satisfy the condition then there is at least one symmetric positive deenite matrix B whose Cholesky factor has the same sparsity pattern as the Cholesky factor of A, for which incomplete Cholesky factorization fails because of a nonpositive pivot.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A Necessary and Sufficient Symbolic Condition for the Existence of Incomplete Cholesky Factorization Xiaoge Wang and Randall Bramley Department of Computer Science Indiana University - Bloomington

This paper presents a su cient condition on sparsity patterns for the existence of the incomplete Cholesky factorization. Given the sparsity pattern P (A) of a matrix A, and a target sparsity pattern P satisfying the condition, incomplete Cholesky factorization successfully completes for all symmetric positive de nite matrices with the same pattern P (A). It is also shown that this condition is...

متن کامل

Effective Preconditioning through Ordering Interleaved with Incomplete Factorization

Consider the solution of a sparse linear system Ax = b when the matrix A is symmetric and positive definite. A typical iterative solver is obtained by using the method of Conjugate Gradients (CG) [15] preconditioned with an incomplete Cholesky (IC) factor L̂ [4]. The latter is an approximation to the (complete) Cholesky factor L, where A = LL . Consequently, the process of computing L̂ relies to ...

متن کامل

CIMGS: An Incomplete Orthogonal FactorizationPreconditioner

A new preconditioner for symmetric positive definite systems is proposed, analyzed, and tested. The preconditioner, compressed incomplete modified Gram–Schmidt (CIMGS), is based on an incomplete orthogonal factorization. CIMGS is robust both theoretically and empirically, existing (in exact arithmetic) for any full rank matrix. Numerically it is more robust than an incomplete Cholesky factoriza...

متن کامل

Locally Optimized MIC(0) Preconditioning of Rannacher-Turek FEM Systems

In this paper Rannacher-Turek non-conforming rotated bilinear finite elements are applied for the numerical solution of second order elliptic boundary value problems. The preconditioned conjugate gradient method is used for the iterative solution of the arising linear algebraic system. A locally optimized construction for an M-matrix approximation of the global stiffness matrix is the first ste...

متن کامل

Incomplete Cholesky Factorization with Sparsity Pattern Modi cation

This paper proposes, analyzes, and numerically tests methods to assure the existence of incomplete Cholesky (IC) factorization preconditioners, based solely on the target sparsity pattern for the triangular factor R. If the sparsity pattern has a simple property (called property C+), then the IC factor exists in exact arithmetic. Two algorithms for modifying the target sparsity pattern to have ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1995